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Abstract. We study the inhomogeneous chiral condensates in the proximity of the chiral tricritical 
point (TCP) of two-flavor QCD. Deriving the Ginzburg-Landau (GL) functional up to the eighth 
order in the order parameter and its spatial derivative, we explore off the TCP and find that critical 
I curves are bent by non-linear effects. In the newly extend GL coupling space, we find the TCP 

being realized as a multicritical point where five independent critical lines meet up. We also present 
general analyses for the energies associated with several higher dimensional crystal structures. 



(N 



(N 



Keywords: QCD, critical point, inhomogeneous phases, chiral symmetry 
PACS: 12.38.Mh,21.65.Qr,25.75.Nq 



43 

Q-i Introduction. - Chiral condensates with possible crystal structures are of particular 

interest in the field of QCD under extreme conditions. Such a partial spatial breaking 
of quark-antiquark condensate may be driven by the quark chemical potential /! serving 
as an external field which favor to produce a net population excess of quarks over 
antiquarks. This was indeed proven to be true in the neighborhood of the chiral tricritical 
point (TCP) in Ref. [1] where the author took the advantage of the Ginzburg-Landau 
^ ! (GL) framework which is capable to make model independent predictions. TCP was 

i/S \ found to be replaced by the Lifshitz point (LP) from which a solitonic chiral phase 

expands. Such phase structure is now known to survive after incorporating the effect of 
current quark mass [2] or Polyakov loop [3]. 

Most studies so far restrict the analyses to the one dimensional crystal structure mainly 
for technical reason, although such state is known to be fragile against thermal fluctu- 
ations [4]. Also moving off the TCP, we might expect other types of inhomogeneous 
states, including the chiral density wave [7], come into play and change the phase struc- 
ture just as in the case of superconductivity under an external magnetic field [5]. 

In this article, we report the extended GL analyses for the chiral crystal phases near 
TCP in the chiral limit based on our recent paper [6]. We aim at the clarification of 
(1) how the inhomogeneous chiral condensate is deformed or taken over by other phases 
once departing from the TCP, (2) whether or not multidimensional crystal structures get 
favored over one-dimensional solitonic state, and finally (3) the structure of the TCP in 
a wider extended GL coupling space. 

Extended Ginzburg-Landau approach. - In order to allow for a minimal description 
of the TCP and possible inhomogeneous states in its neighborhood, we need to expand 
the GL functional up to sixth order in the order parameter and its spatial derivative. The 
expression was found as [1]: 



CO = ^M(x) 2 + f (M(x) 4 + (VM) 2 ) + ^ (M(x) 6 + 5M 2 (VM) 2 + ^(AM) 2 ) , (1) 
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FIGURE 1. (Left panel): The GL phase diagram for Eq. (1), computed at the mean-field level with 
restricted to one-dimensional structures. (Right panel): The phase diagram after including the effect of 
eighth order terms, 8(0. For comparison, the critical lines enclosing the inhomogeneous phase in the 
absence of 8(0 (those in the left figure) are repeated by thin dashed lines. 



where M(x) ~ (qq(x)) is the chiral condensate, and OC2, OC4 and Ofe are the relevant GL 
couplings. If CC6 > 0, We can draw a phase diagram in the two dimensional GL parameter 
space ((X2, OC4). It can be shown that the thermodynamics only depends on the sign of OC4 
and dimensionless ratio / '((^[oc^ 1 ]) where (^[o^ 1 ] means OLiO^, i.e., the value of 0C2 
measured in the unit of CC 6 l . Due to this scaling property, all the critical lines in the 
(fl^Gg ]> sgn(a4)a|) -plane should be expressed as straight lines, and any two of them 
never intersect one another once departing from the origin. The phase diagram for this 
case is shown in the left panel of Figure 1 . We find a Lifshitz TCP at the origin; the 
region labeled by M(z) is the phase of inhomogeneous chiral condensate characterized 
by Jacobi's elliptic function M(z) = y/kvsn(kz, v) with v being the elliptic modulus. 

Exploration off the TCP. - In order to explore a new phase structure which may show 
up away from the TCP within the GL framework, we need to go beyond the minimal GL 
described above. The term which should be added at eighth order is [6] 

8ca = ^ (M 8 + 14M 4 (VM) 2 -i(VM) 4 + fMAM(VM) 2 + ^M 2 (AM) 2 + i(VAM) 2 ) . 

(2) 

Coefficients appearing in front of each term are all proportional to otg, and can be ex- 
tracted by evaluating the quark loop diagrams. 0^ should be positive for thermodynamic 
stability, while can be either positive or negative. We first focus on the case ctg > 
which has close connection to the analysis at the sixth order, and later come back to the 
case «6 < 0. With inclusion of this 800, we can see that the thermodynamics now is a 
function of two independent parameters, 0C2 and OC4. In the right panel of Figure 1, the 
phase diagram is depicted in the ((X2, sgn(a4)a|) -plane. In this case two critical lines 
surrounding the phase are bent by the nonlinear effect from the eight order terms. In fact 
it is possible to derive the following expansions for the formation of single soliton (the 
dashed curve on the left side) and for the melting of condensate (the one on the right 




side, red online): 

cd = ^ocf (1 + ffa 8 a4 + H (as^) 2 + 6 ((a 8 «4) 3 )) , 

(3) 

a 2 = |a| (1 + ^a 8 «4 + jm^®*) 2 + & (("s^) 3 )) • 

Here «2 and ctg in the above formulas mean o^fo^ 1 ] and as[a|] respectively. However, 
apart from these nonlinear bending of critical lines, no new phase structure appears. 

Multidimensional crystal structures? - We have limited the above analyses to the 
condensate having its spatial modulation only in one direction, i.e., one-dimensional 
modulation. With this restriction, it is possible to solve the Euler-Lagrange equation 
analytically leading to results demonstrated in Figure 1. Here we lift the limit, and 
explore the possibility of multidimensional crystal structures using a variational method. 
The states we consider are the extensions of one-dimensional Larkin Ovchinnikov (LO)- 
type condensate (Mid-lo(x) = V2Mosm(kz)) to higher dimensions, 

M 2D -lo(x) =M (sin(fct) +sin(£y)), 

ft (4) 
^3D-lo(x) = y §M (sin(fac) + sin(ky) +sin(fe)), 

where k and Mq are to be optimized for each state. Figure 2 shows the comparison of the 
spatial average of order parameter (M(x)) (Left panel) and the averaged potential (d)(x)) 
(Right panel) for the several states along the section = — 1 in the left panel of Figure 1. 
"FF" refers to the Fulde-Ferrell (FF)-like state characterized by Mid-ff( x ) = M$e z in 
which the time-reversal symmetry as well as rotational symmetry is broken. We see that 
the magnitude of condensate is decreasing function of dimensionality and the energy 
density is accordingly in the opposite tendency. Similar results for two dimensional 
condensates at T = were recently obtained in the NJL model [8, 9]. The inclusion of 
eighth order terms in 8(0 does not change this situation so that one-dimensional soli tonic 
state is most favorable for all values of a^a^. 

We can derive the analytic formulas for the dependence of free energy on the number 
of spatial directions in which the condensate varies by expanding the GL functional in 
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FIGURE 3. The GL phase diagram for OLf, < computed from Eq. (1). Dashed lines express the second 
order phase boundaries, while full lines do the first order ones. Dot-dashed line in the inhomogeneous 
phase represents the first order phase transition which would have been realized if the inhomogeneous 
phase was not taken into consideration. 



the vicinity of condensate-melting line. We first set the higher (d-)dimensional analog 
of the LO and FF states as 

M dD _ FF (x) = ^M (e tol + e fa2 + • • ■ + e lkx "), 

r (5) 
^dD-Lo(x) = ^M (sm(kxi) + sm(kx 2 ) ^ \-sin(kx d )). 

Then at the proximity of the line of vanishing condensate, we arrive at the following 
analytic formulas for (o(x)) via expanding in Mq after optimizing the wavevector k 
over Mo: 

(«dD-Lo(x)> = (f 

(C0dD-FF(x)) = (f 

From these we see that, in the case of LO-type condensates the energy density is in the 
order of dimensionality of spatial modulation. In the FF case, the difference becomes 
visible only at sextic order, and the dependence on dimensionality is not monotonic; 
the 2D structure is most favorable within this class of condensates. From these analyses 
we find that increasing the directions of spatial modulation is accompanied by a large 
kinematical energy cost, leading to the conclusion that multidimensional crystals are 
unlikely to be realized in the neighborhood of the TCP. 

The structure of TCP. - We now explore a new regime Ofe < 0- The corresponding phase 
diagram is depicted in Figure 3. We see that the Lifshitz tricritical point for Ofg > splits 
into four multicritical points; the critical point (CP), the Triple point (T), the Lifshitz 
bicritical point (L), and the critical endpoint (E). Their (0C2, OC4) -coordinates are listed 
in TABLE 1. All these points scale as (# | ocg | 3 /ag,#a^/as) so that they are all going to 
degenerate into the origin as (X(, — > ~, smoothly continuing to the Lifshitz TCP in the 
positive «6 side- Among these critical points the most intriguing is T where three first 
order phase transitions meet. At this point, three different forms of broken-symmetry 
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TABLE 1. The coordinates of locations of the multicritical points. 



Label 


Type 


Coordinate (054 


[« 6 3 /«8 2 ]>02 [a 6 2 /«s]) 


CP 


Critical point 


(-1/27,1/3) 


(analytical) 


T 


Triple point 


(-0.016,0.27) 


(numerical) 


L 


Lifshitz bicritical point 


(0,5/18) 


(analytical) 


E 


Critical endpoint 


(0.014,0.16) 


(numerical) 



phase compete and coexist. These are; the state with homogeneous chiral condensate, 
the state with smaller chiral condensate denoted by "#SB2", and the inhomogeneous 
chiral condensate labeled by M(z) . It is possible to extract numerically the dimensionless 
ratios of physical quantities that are to be realized when the point T is approached: 

(cfccHT-tMOtSB)' M( X SB) 'yfJMWjl 

where £p is the modulation period. These ratios characterize the triple point (T) and they 
are universal and model-independent in the sense that they do not depend on any of GL 
parameters including the magnitude of | OCg | , and a&. 

Summary. - In conclusion, we explored off the TCP based on the GL functional ex- 
panded up to eighth order in chiral order parameter and its spatial derivative. As a con- 
sequence, we found that the critical lines surrounding the inhomogeneous phase are bent 
because of the nonlinear effect coming from eight order terms although the qualitative 
phase structure remains unaffected. We also examined the possibility of higher dimen- 
sional crystal phase, but it turned out that the kinematical cost to make the condensate 
modulated in several directions is so high that it is unlikely to have multidimensional 
chiral crystals near the TCP. Finally we extended our analyses to a new regime < 0. 
We found that the TCP splits into four multicritical points. It would be interesting to 
search for a thermodynamic physical quantity in QCD to which 0^, is sensitive. The 
isospin chemical potential or the chiral chemical potential might serve as one of such 
sources. The investigations in these directions deserve further efforts in future [10]. 

Acknowledgments. - H. A. thanks the organizers of QCD@Work2012 for giving him an 
opportunity to give a talk. The authors thank D. Ishibashi for the fruitful collaboration. 



REFERENCES 

1. D. Nickel, Phys. Rev. Lett. 103, 072301 (2009) [arXiv:0902.1778 [hep-ph]]. 

2. D. Nickel, Phys. Rev. D 80, 074025 (2009) [arXiv:0906.5295 [hep-ph]]. 

3. S. Carignano, D. Nickel and M. Buballa, Phys. Rev. D 82, 054009 (2010) [arXiv:1007.1397 [hep- 
ph]]. 

4. G. Baym, B. L. Friman and G. Grinstein, Nucl. Phys. B 210, 193 (1982). 

5. S. Matsuo, S. Higashitani, Y. Nagato, and K. Nagai, J. Phys. Soc. Jpn. 67, 280 (1998). 

6. H. Abuki, D. Ishibashi and K. Suzuki, Phys. Rev. D 85, 074002 (2012) [arXiv:1109.1615 [hep-ph]]. 

7. E. Nakano and T. Tatsumi, Phys. Rev. D 71, 1 14006 (2005) [hep-ph/041 1350]. 

8. S. Carignano and M. Buballa, arXiv:l 1 1 1.4400 [hep-ph]. 

9. S. Carignano and M. Buballa, arXiv: 1203.5343 [hep-ph]. 

10. Y. Iwata, H. Abuki and K. Suzuki, arXiv: 1206.2870 [hep-ph]. 



